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SUMMARY 

An analytical method is presented which approxi¬ 
mates the flight-path deviation of a variable-mass, 
variable-inertia , spin-stab dazed vehicle under the 
influence of ( a) initial pitching motions, (b) thrust 
ro isalinements, and (c) principal-axis misalinemeni , 
or dynamic imbalance. I his method includes a 
first-order analysis of the effect of jet damping on 
the resulting flight path. 

An extensive analysis was carried out on an IBM 
70j electronic data processing machine , by integrat¬ 
ing numerically the six-degree-of-freedom equations 
of motion of a thrusting body with variable mass and 
inertia. This computer analysis was intended to act 
as a check on the method developed herein, and also 
to present other investigators with exact numerical 
answers with which to compare their own results. 
The results from the analytical solution are com¬ 
pared with the results from the numerical integration f 
and the results are found to agree well in all cases 
where no large change in inertia ratio was involved. 

The results are presented in terms of the angular 
deviation in space of the velocity vector from its 
undisturbed orientation , as a function of spin rate, 
for each of the disturbing parameters. It is belie red 
that this velocity-vector deviation as a function of 
spm rate is a much more intuitive indication of the 
dispersion than the body attitude in space, which is 
the m ost fre que fitly used pa ram e te r in sp i n -sta b Hi za- 
tion discussions. 

INTRODUCTION 

With the fairly recent advent of very-high- 
altitude rocket vehicles, the engineer has become 
facial with the difficult problem of adequately 


stabilizing his vehicle in the region where the 
atmosphere is so thin that conventional aerody¬ 
namic controls are rendered useless. One possible 
means of stabilization is the use of reaction con¬ 
trols in place of the aerodynamic controls. How¬ 
ever, in view of the diminutive size of some of the 
present-day probes and satellite injection stages, 
a control system of this type would impose a 
severe weight penalty on this final stage, lienee, 
some other means of pitch and yaw attitude 
stabilization must be considered. 

Ballisticians have long known that a spinning 
shell deviates from a predetermined flight path 
much less than a nonspinning she'll. It thus 
seems natural to attempt to attitude stabilize a 
rocket vehicle by spinning it and, indeed, this 
method has been successfully used already in the 
Vanguard, Kxplorer, and Pioneer vehicles, and on 
Scout and other vehicles currently being designed. 
Spin, however, is not a cure-all: spinning a vehicle 
gives it a constant attitude in inertial space, but 
not stability along the constantly changing tan¬ 
gent to tiie flight path. Hence, if a large change 
in flight-path angle is involved, the method of 
spin stabilization may not be applicable. Spin 
stabilization is very effective, however, in the case 
of a satellite launcher where the last stage, the 
injection stage, is final at nearly horizontal flight- 
path and attitude angles and it is desired to burn 
out this stage in such a manner that the final 
velocity vector will be nearly horizontal. 

Some investigators in this field (for example, 
refs. 1 and 2) have chosen as a criterion for spin 
stabilization the attitude of the body, defined by 
tli(‘ two Puler angles in pitch and yaw. It- is fell, 
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however, that a knowledge of the vehicle attitude 
alone is not the best criterion for determining the 
effects of various extraneous forces, because the 
attitude alone does not indicate the final orien¬ 
tation of the velocity vector at stage burnout. 
This velocity-vector orientation is believed to be 
a better indication of the dispersion effects sought. 
Consequently, the data presented herein are given 
as deviation of the velocity vector from its undis¬ 
turbed orientation as a function of spin rate* for 
several external disturbing forces. 

The two methods in most general ust‘ for deter¬ 
mining dispersion ot a spinning vehicle, with the 
body attitude as a criterion, are the methods of 
Nicola ides (ref. 1) and Jarmolow (ref. 2). Xico- 
laides determines the dispersion in terms ot angle 
of attack a and angle of sideslip (J -his parameter 
is essentially % a 2 +>\ He also includes aerody¬ 
namic effects, which are neglected hero. However, 
he considers the vehicle mass and moments of 
inertia, as well as the applied moment and forward 
velocity, to be constants. His results, therefore, 
cannot be used to study the dispersion of a spin¬ 
ning rocket during its thrusting phase. 

darmolow, on tin* other hand, considers the case 
of variable mass, moments ot inertia, and applied 
moment ; but lie diseusses only the 1 body attitude. 
Also, analytical results are given only for the. pitch 
and vaw rates in a body-axis system. Numeiical 
integration must still hi* list'd to find the pitch and 
yaw Euler angles. 

The method presented herein utilizes exponen¬ 
tial approximations for the variation ot the ratio 
of the applied moment- to 1 he pitch moment of 
inertia, and for the ratio of thrust to linear mo¬ 
mentum, which appears as a parameter in the 
translation equations. 1 hese approximations 
make possible analytical solutions for pilch and 
yaw Euler angles, and also expressions defining the 
angular deviations of the velocity \cctor. I hese 
expressions include linear damping in pitch and 
yaw. The following additional assumptions are 
’made: tin* ratio of roll inertia to pilch inertia, is 
constant: the vehicle is rigid and lias rotational 
svmmetrv; then* is no roll moment; only small 
angles are considered: the motion takes place in a 
vacuum. This study was initiated to determine 
die effect of various disturbing parameters on the 
flight path of a solid-propellant rocket motor 
similar to the final stage* of the Seoul. 


SYMBOLS 


constants 

c ' pitch and yaw damping coefficient, 

ft-lb-sec 


r 



(\,<y i, 

(\ x constants 

H„ 

angular-momentum vector of ve¬ 

H, 

hicle with respect to center of 
gravity in coordinate 

system, slug-fV/sec 
angular-momentum vector of ve¬ 

II, 

hicle with respect to center of 
gravity in x n y n z c coordinate sys¬ 
tem, slug-f t “/sec 

angular-momentum vector of ve¬ 

Li , J h j k (, 

hicle and jet with respect to 
center of gravity in x rf y n z r coor¬ 
dinate system, slug-ft-/See 
unit vectors along the AV, 3 

LJr.k 

and Z f -axis, respectively 
unit vectors along the AV, )V, 

/.vj / vi y. 

and Z r -axis, respectively 
vehicle moments of inertia about 

I 

AV, >V, and Z„-axis, respec¬ 
tively , slug-ft - 

J y or // (where Iy—Iyj, slug-ft J 

I } Av.?! 

,/ vv vehicle products of inertia about 

j ^ 1 j ^ 2 

AV, V 6 -, and Z,-axis, respec- 
tively, slug-ft- 
constants 

/ 

length from vehicle center of grav¬ 

M 

ity to nozzle exit plane*, ft 
disturbing-moment vector, d/.vL 

M 

Myk n ft-lb 

magnitude of disturbing moment, 

in 

ft-lb 

mass of vehicle, slugs 

in,, 

element of mass, slugs 

u J ! 

-*• 1 * U 

at. t — U, 1/sec- 

S 1 i’¬ 
ll 

„ radians/scc 

/v/z 

angular-velocity components of 

T e 

vehicle about AV, IV, and Z f - 
axis, respe(*tively, radians see 
unless otherwise indicated 
distance from center of gravity to 

r e 

jet exit, r<. x \ < .+ r, iY ]c+ f’r./Mc, E 

magnit ude of r,, Cl 
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t 

time, sec 

T 

t hrust, lb 

V 

vehicle velocity vector, ft/sec 

V 

magnit ude of V, ft /sec 

V 

jet exhaust velocity, bV.\-i r + 
Ve,rJr+Vr.zK, ft 

r. 

magnitude of V 0 ft/sec 

J h 1 }j hy~ t> 

coordinates measured along body 
AV, ) V, and Z,-axis, respec- 
t ively 

A /l, ) hy h 

body-fixed system of orthogonal 
axes, origin at body center of 
gravity 

j c 1 // Cy~r 

coordinates measured along AV, 
>V, and Z ( .-axis, respectively 

A C, ) 

system of orthogonal axes which 
pitches and yaws with vehicle 
but does not roll, origin at body 
center of gravity 

3 u Hi a 

coordinates of (‘lenient of mass 
causing dynamic unbalance 

Cl 

angle of attack, radians 

V 

angle of sideslip, radians 

y 

orientation of velocity vector in 
vertical plane, radians 

6 

thrust misalinement angle, radians 
(except as noted) 

e 

principal-axis misalinement angle, 
radians (except as noted) 

A 

orientation of velocity vector in 
horizontal plain*, radians 

(7 

total orient at ion of velocity vector, 
radians (except as noted) 

0 

pitcli Euler angle, radians (except 
as noted) 

* 

yaw Euler angle*, radians (except 
as noted) 


angular velocity of nonrolling sys¬ 
tem in inertial space, a> rr j ( .4- 
o>zM n radians/sec 

(*>' lV) W Zc 

Subscripts: 

components of os r along )V and 
Z r -axis, respectively, radians/sec 

0 

denotes quantity at zero time 

X.VyZ 

components parallel to AV, 3V, 
and Z A -axis, respectively 

AV,3 V,Zc 

components parallel to AV, JV, 
and Z r -axis, respectively 

A dot over 

a quantity denotes differentiation 

with respect to 

t ime. 


ANALYSIS 

The mathematical analysis of a spinning rocket 
differs in several respects from that of a spinning 
projectile. Jn the first place, the rocket has mass 
and inertia properties which vary with time, 
whereas those of a spinning projectile are constant. 
Secondly, the rocket can have a moment due to a 
thrust misalinement. This moment, which is 
fixed to the spinning vehicle, gives rise* to motions 
not found in the analysis of projectile motion. 
The motions due to an initial pitcli or yaw rate 
and a dynamic unbalance 1 an* common to both 
types of bodies, but the analysis becomes more 
involved when a variable-property vehicle is dealt 
with. 

The (‘((nations of motion programed on the 
IBM 704 electronic data processing machine for 
comparison purposes in t his st ud v are t he standard 
six-degree-of-freedom ('((nations found in most 
aeronautical textbooks (for example, ref. V. 
The ('((nations programed, along with the various 
transformations that were used, are presented for 
reference in the appendix. 

An approximate analytical method of predicting 
the deviation of the velocity vector due to the 
effects of thrust misalinement, initial pitch rate, 
and principal-axis misalinement will be developed 
in this section. Two coordinate systems will be 
used in this analysis. The first is a 

body-fixed coordinate system with axes fixed at 
the center of gravity of tin* vehicle and the 
coordinate x tl extending along t he body longit udinal 
axis. The second coordinate system (./>,// 0 ~r)? hi 
which the equations of motion will be derived, is 
referred to a body-fixed system of axes which can 
pitch and yaw with the vehicle but does not roll 
wit hit. The t wo sets of coordinate axes are alined 
at /. — (). (See fig. 1 .) The pitch Euler angle 0 
and the yaw Euler angle ^ are assumed to lx* small 
so that the older of the Euler rotation is immate¬ 
rial. This assumption of small angles simplifies 
tin* integration for 7 and A, the flight-path angles 
in the vertical and horizontal planes, respectively, 
and also permits these flight-path angles to 
combine vectorially to give <r, the angle between 
the velocity vector and some fixed reference line, 
taken hen* as the longitudinal axis of the vehicle 
at zero time. 

The angular-momentum vector of the vehicle 
in the f b ,y h ,z h system instantaneously fixed in 
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Fk.i kk 1 -Coordinate system and definition of symbols used in analytical development. 


inertial spare can be written 

Hj, ~ - [T\ p l\ y(j Ixz f ! h> ~r {fy ( J 1 Yz f lx yP)\ b 

f\'/.P I ( 1 ) 


Now, in the non rolling system (x n y cf z c ) f which 
has angular velocity in inertial space, 

C0,= (Q,Wy Cf O>z c ) = U)Ycic-\ UzMc (2) 

The conation of motion about (he vehicle center 
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of gravity is 



(:*) 


The rate of change of angular moment urn with 
respect to the vehicle center of gravity including 
jet effect is 


(IHr 

' (It 


(i n r 

(It 


ni r e X ( 




V,) 


(4) 


where wXH r takes into account the rotation 
of the* coordinate* systeun. 

it is usually conventional in missile* work to 
choose* the body axe*s as principal axes and to 
assume a symmetrical vehicle, so that 


I y — Iy— I 
IX 1 ■ — I yy — I X y ~ 0 



With the *se* assumptions, expiation (1) may be* 
written in the* nonrolling system as 

H c r ~ IxpK + Iu>Yrl c 4 ItozMc (ti j 

If the* further assumptions are* made that tlie 
jet is located on the AVaxis and the* je*t <*xit 
velocity redative to the* nozzle <*xit is dire*(*te*el 
along the* negative* AVaxis, tlien 


and 


’ (. X c * t 

l r Yc~~ 1 r.zc - o 




C. X c i r 

ff.Yc ? c, Y.t b 



Expansion of equation (3), with the aid of sub¬ 
stitutions from eepiations (2), (4), and (0), yiedds 


(/ v P “t - / xp )b A (/uj Yc -r 7a > Vc -f- 1 x ])u j, r -p M f iw > V j r 

— (I(C Zc /<o, c — I yp to y c + th fi 00 , f ) k c — M (9) 

if damping in roll is negl<*et(*<i (that is, if f x - 0), 
equation ( 9 ) becomes, in component form, 

fxi> = Mx (10a) 

7to rc ^/ v />to, r + c'co rc - - \f y (10b) 

l<J?Zc IxP&Yc~t 6 Cjd, c 4/, ( 10c) 


whe 


re 


c f = mrl I 


The quantity c'/I( = c) will be considered constant 
for this analysis. 

If the* assumptions are made that tho disturbing 
moment lie*s in the vertical plane* at ze*ro time 
and that the*re* is no moment about the* roll axis, 
t hen 



-l/.v = 0 'j 



My— M cos pt >- 

(11) 


Mx — M sin pi J 


it ions 

(10) become 



Zv/b-O 

(1 — si) 

/w Y c 

4 ~ fxpuzc-he'uye — Af cos pt 

(121.) 


— F\pO)y c y (Af Sill J)t 

(12c) 


Eejuation (12a) immediately show's that 


Ia*t 


and 


p -( V)nst ant jt 0 


[x 

I 


At 

I 


Pv=H 


My 


k! 


3 ) 


14) 


It is assumed that the ratio f x jl is constant: 4/ y 
and k are* also constants, found by fitting an 
exponential to the* function Mjf as shown in 
(equation (14). Then e*quatiems (12b) and (12c) 
be*c*ome* 

a ) Yc J r noi Zc -rC^ Yc ^ Af 0 e k( cos pj (loa) 


co,,— nw Yc -\-cu Xc ^ MM 1 sin pj 


(lob) 


'These* equations earn be integrated simply, and 


with the* 

initial conditions co rr (())-- 

Yc (^ 7 .i 

(0)=0 

give* 




Uy c — («r 

— fV)e _c ' e*os nt — (\c~ a 

sin nt 



~r-(y kt COS pJ-r(\c H 

' sin pj 

(1 (ni) 

&Zc — V'Ci 

, t 0 — ( V)r ~ cl sin nt- (\< ~ ct 

ceis nt 



- l - ( y * sin pj — (\e 

kt cos pj 

(101.) 

where* 





M 0 (p„ — n) 

(c + k )“ 4 - {p ti i 

y 

( 17 a) 
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.\fjc\-k 


(c -4- k) ~ -f- {p<> n ) “ 


(17b) 


Now, by virtue of the original assumption that 
the pitch and yaw motions are small, and as a 
consequence of the coordinate system used, the 
following expressions can be written: 

6— J'tPypif (ISa) 

i J'uzftii (hsh) 


Substitution of equations (16a) and (16b) into 
these integrals with the initial conditions 
6 ( 0 )=\p (()) — 0 gives 

S—A(e kt cos pJ—\)-\-Iie kt sin pj 

— [){e~ ct cos ///— 1 ) -\-Ke~ ct sin nt (19a) 

\p~ Ae kt sin pj—li{e kt cos pj— 1 ) 

— })(~ c! sin nt — FAe~ ct (‘os nt— 1 ) (19b) 


where 


C 1 k-<\l> 
F \ Pi 


(20a) 


If— 


F-\ Pi 



n 


r(a'}v, — ( 7) —*. 
c~-\- tr 


(20c) 


considerations. The equations defining the rota¬ 
tional motion of the vehicle are completely inde¬ 
pendent of the translational motion. (See, for 
example, eqs. (Al) to (AG).) Hence, for given 
inertia characteristics, spin rate, and disturbance, 
the rotational motion of the vehicle is completely 
defined. Now, in order to get the actual transla¬ 
tion in space, the body attitude parameters must 
be combined with the translation parameters; 
for example, equations relating quantities such as 
thrust and mass ratio are integrated to give the 
translational motion. As an illustration, it is 
certainly possible to have a spinning vehicle point¬ 
ing vertically upward and moving horizontally at 
ignition. Then as the motor tires the velocitv 
vector climbs above the horizontal. Thus, the 
final direction of the velocity vector has changed— 
the dispersion as defined here is reduced—but the 
body attitude has remained constant. An effect 
such as this one is not evident when only the rota¬ 
tional motion is defined. 

Jn order to determine y, tin 4 orientation of the 
velocity vector in the vertical plane, reference is 
made to the normal equation of motion of figure 
2(a). Since a —6 — y, if only small angles are 
taken ii to account this equation becomes 


fly, T T _ , T . 

,. + ■ ,-.7=—^ 04- .rOCOApJ 
(I t m \ m l in l 


( 21 ) 


In a sii filar manner a differential equation for X, 


K 


1 c( 


»2 _ L, 
[ 


/c 



These equations give time histories of the body 
attitude angles in an inertial coordinate system 
taken to he coincident with the body coordinate* 
system at zero time. These angles 0 and arc* 
the angles eommonlv used to discuss the effects of 
spin on the stability or dispersion of a vehicle. 
However, body attitude* alone does not adequately 
describe tin* complete motion in space. The act ual 
dispersion of the vehicle, defined herein as the 
angle* between the velocitv vector at zero time and 
the velocity vector at burnout, is a much more 
meaningful parameter because* it depends on the* 
translation of the* vehicle, and lienee on the quan¬ 
tities which define its translation history, as we*ll 
as on the rotational parameters. The* pm‘c*eling 
argument can be substantiated from physical 



mVy - Sin 8 cos p Q t COS a + f COS 8 sin a 
mVy = T 8 cos p Q t + 7"(0-y) 

(a) Vertical plane. 

ITco rk 2. Equal ions of mot .ion. 
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' orientation of the velocity vector in the hori¬ 
zontal plane, can be obtained from figure 2 (b): 


ti \. r x r , , r . 4 

<lt + ml' x ~ "mV* 1 mr 5sll| yC 


( 22 ) 


Since only the changes in 7 due to external dis¬ 
turbances rather than (lie natural change in flight 
path due to gravity are of interest, the effect of 
gravity is omit ted. 

Now, in order to integrate equations ( 21 ) and 
( 22 ), the quantity TjmV must bo put into a form 
that is easilv integrated. It can be shown that a 
plot of T/mV as a function of time can be approxi¬ 
mated to a sufficient degree of accuracy by an 
exponential function. Thus, define. 


T 

/at 


T =kyt k 2t 


(2:ft 


and 


r ki 


k 


k. 


In tlie process of integrating equations ( 21 ) and 

( 22 ), tlie quantity c* 2 V A ‘ e * 2t appears in several in- 
staiu*es. Therefore, to further simplify Ihe inte¬ 


gration, it was assumed that. 




<\e r * 1 





i* T cos 8 


m V X = T sin 8 sin p 0 1 cos fi + T cos 8 sin 
m V A - T B sin p Q t + Tj ,S 

(h) Horizontal plane. 

Fkjukk 2 . Concluded. 

If equation ( 2 d) is substituted into equations 
( 21 ) and ( 22 ) along with tlie expressions for 6 and 
\p from equations (It)), integration of the resulting 
linear differential equations, utilizing equal ion (24) 


and the initial conditions 7 ( 0 ) A( 0 ) 0 , yields 


7-r 


ke 



Ae u ' 21 '((',+/•) -/>„lie “■'*' * 1 ' 8e v *‘< 

(TiA-ky+pl + r>+ P i_ 

cos pj + 

Ap,,e' 

^1 

•so Ly 

cNo; ? 

i_i 

ri)((J.,— c)e l< ' 2 ~ cU -\-p’ne^~' 2 ~ cu ~ 

sm p u i 1 ' • _ - 

L (O-OM-m- J 

cos nt —T 

HI 


+ /^ ( ^ |tW (f' 3 - 1 /*) 
2 + k) 1 +]>l 


{('t-cy+n* 


» 

si 


sin nt 


VAf* D<f*n r A(C 2 +le)-Bp 0 I)((' 2 -c)+En A 1) d(' 2 , 

L (', T\ J L ( (T 2 +ky+i>i (r 2 -c o 2 +w* V 2 c 2 cj+rijj 


- k.,t _ 

\=er ke 




— p„Ae «*+*>«_ Ii«’ 2 +k)e t( ~ 2+ku Po 8e c *' 


«v\ -ky+pt 


p„8e c ' 2 ' I . f* 

-tswJ+L 


Ae^+v'iC.+k) - pJie iC * + * H 

(71+1)* I pi 


+ 


<w*~\ . ,,rr 

_■ - sin pj I _ - COS III — 

CM piA L (O-r) 2 I «* J L 


CU (C.,- c)+nEe u v r " 
(O-r)M «* 


si 


sin nt 


rji/v E£*n_r-p^\-i{((’^-k) _ _p„8 ni)-{<? a -c)E n k i \ ( 

L r, r, J L (.r 2 +ky~\ pi r 2 \ pi (A-c)M « 2 r 2 Vj/ 


021963 -02 -2 



8 


TECHNICAL REPORT R-110-NATIONAL AERONAI* 1'ICS AND SPACE ADMINISTRATION 


Tii (Ml, because* of the* assumption of small angles, 
7 and X ran be (*oml)incd vectorially to give* a, the 
angle between the fixed reference line and I lie 
velocity vector: 

<r-y?+X 2 (27) 

A disturbing parameter consisting of a princi¬ 
pal-axis misalincmcnl can lx* handled in tin* fol¬ 
lowing manner. Assume that the mass element 
causing the dynamic unbalance 1 can be placed oil 
the center line of the vehicle a distance >/ u from 
the center of rotation, and that it lias a mass 
in u , as in sketch l. The moment about the renter 



Sketch 1. 


of gravity due to this mass element, is 

A/. pU„y«Mu (2N) 

Now, tin* produc t of inertia due to this mass is 

/.vi* - (29) 

Hence, the moment can lx* written in terms of tin* 
product of inertia: 

u-pHxY m 

It may lx* desirable to express the product of 
inertia in terms of the principal-axis misaline- 
menl angle e. From the properties of moments 
of inertia it is known that the* angle e between 
the principal axis and the axis about which the 
moments and products of inertia are computed 
is given by 

lau'ie— vv , (HI) 

* x / 

If, as will be the usual cast*, the product-of-in(*rtia 
term is due to construction technique or placing 
of instrumentation, the product of inertia will be 
constant, the small change* due to tlx* center-of- 
gravity shift during burning being neglected. In 
this case, e will lx* a function of time*. The time* 
history of e can bo computed from equation (31) 
since* / v and / are known at any time. Then 


with the assumption c<< 1 , substitution of equa¬ 
tion (31) into equation (30) give*s 

'/ /*'/(i -f) m 

If this expression is used in expiation (14), and M a 
and k arc* found by a curve* lit , solution of expia¬ 
tions (25) and (20) will yield histories of 7 and X 
due to a principal-axis misulincment e. 

It should lx* noted that, to the degree of approxi¬ 
mation usexl herein, the* (‘fleets of the various 
disturbing parameters can lx* combined linearly 
to yield the total dispersion. For example, equa¬ 
tions (25) and (20) may he* written as a linear 
combination of terms each of which is a function 
of ojyc.i, M 0 , Dr 8 but not of any t wo of the variable's. 

It should lx* further noted that, again to llu* 
degree of approximation assumed here, the angles 
7 and X, and consequently the total dispersion a, 
vary linearly with a ) YCt0 and Aif taken sepa¬ 
rately, so that it is not necessary to make* a 
separate calculation for each value* of moment or 
pitch ^ate desired. This linearity is supported 
by tlx calculations performed on the* IBM 704 
romputcT, as illustrated in several of the* figure's. 

RESULTS AND DISCUSSION 

A comparison of the results of the* numerical 
integration (on tlx* IBM 704 computer) of tlx* 
exact -*epiations of motion given in fix* appendix 
with tlie results of the analytical solutions de- 
veloped herein is made*- in figure's 3 to 10 . The* 
qimntuie's 0, \p, 7 , X, and a we*re calculated from 
('(Illations 09a), (19b), (25), (20), and (27), 

respectively. The* constants used are given in 
tables I and II. Figures 7 to 10 are plotted for 
/ = 30 seconds (burnout). These figure's serve a 
double* purpose. Since not everyone has access 
to an IBM 704 or to a computer of sufficient 
capability and storage capacity to permit the 
progra ning of a problem as complex as a six- 
degree-of-freedom rigid-body analysis, many in¬ 
vest igi tors have probably developed analytical 
results without being able* to clieck the* accuracy of 
tlie'ir 1 let bods. Thus, the figure's presented can 
seTve to ehe*e*k the'se* me*thods and indicate regions 
of validity and applicability. At the same time, 
the figures serve to illustrate the* accuracy of the 
analytical method developed herein. As <*nn be* 
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TABLE I. PHYSICAL AND PERFORMANCE 
CHARACTERISTICS OF ROCKET MOTOR 
CSKI) IN CALCULATIONS 


Roll moment of inert ia, slug-ft 2 . 5.25 0.18/ 

Pitch (yaw) moment of inertia, slug-ft 2 _ ILL 15 0.88/ 

Mass, slugs _ ..... . _ 18.05 0.18/ 

Applied moment, ft-lb (5 in radians) 10.7105 ; 1505/ 

Thrust, lb ... . 8,000 

Burning t ime, see . . . 80 

Init ial velocity* ft/see .. . 10,000 

Final velocity, ft/sec _ 10,810 

(mm] mted const lints (for 0.10° t li rust misalinement): 

hi.. _. 0.01588 

k>_ . 0.0200 

I. _ ... 0.008 

(\ ___ ....... 1.880 

r 2 _ . . ... o.o5o 

5, radians.. 0.001745 

co» fM radians, sec .. .... ...... 0 

lx!I . -- --. ..... 0.104 

4C, 1 sec-.. .. . . .. 0.5470 

k ........... .. 0.0502 


TABLE II. PHYSICAL AND PERFORMANCE 
CHARACTERISTICS OF SPHERICAL R( )CKET 

MOTOR 


Roll moment of inertia, slug-ft 2 
Piteh (yaw) moment of inertia, slug-ft 2 
Mass, slugs . 

Thrust, 11) 

Burning t inn*, see 
Moment, ft-lb 
Initial velocity, ft /sec 
Final velocity, ft,'see 
(Computed constants: 

4*i - - - - -. 

k> _ __ . . 


888.0 

408.0 

10.88 


0 10 . 0 / 
0 10 . 0 / 
8 0.51/ 

8.400 
80 

5 | 0.01/ 
10,000 
20,050 

0.02125 

0.0700 


4*. __ 


0.200 


C,_ 1.280 

C 2 ___ _ .. .... 0.107 

5, radians _.. . . 0 

ojyc . o , radians/sec ...... 0 

Ix/I- ___ ____0.08 

M„, 1/sec 2 __ _. . . _ 0.0124 

k ... .. . 0.042 


soon from a general perusal of figures 8 to 15. tho 
accuracy obtained is quite good in the range 
invest igated. 

Figure 3 shows tho pitoli Euler angle 0 plotted 
against the yaw Euler angle for the last few 
seconds of burning of the rocket motor, under the 
action of a moment produced by a thrust misalinc- 
ment 5 of 0.10°. This figure' shows, essentially, 
the attitude of the body iu space during this time. 
It can be seen that, although the analytical method 
predicts frequencies of nutation and precession 
which are slightly high, the amplitudes of both 


0,deg 


^,deg 


Finnic 8. Undamped body altitude'. />„—80 rpm; 

V--0.KI 0 . 



motions, which are tin* parameters of primary 
interest, are predicted rather well. This curve 
does not include the effect, of jet damping. 

Figure 4 shows time histories of the deviation 
of the velocity vector a. Tho analytical method 
predicts results which are slightly conservative. 

In figures 5 and b are plot toil time histories of 
the pitch Euler angle 6 and tho yaw Euler angle* \p, 
respectively. Hero it will be noted that, while 



0 4 8 12 16 20 24 28 


Time, sec 

Florae 4. -Deviation of velocity vector as a function of 
time, for a thrust misalinement angle 6 of 0.10°. Un¬ 
damped body; p 0 — 80 rpm. 
























Deviation of velocity vector, <r, de 





Deviation of velocity vector, cr, deg 
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Figure 8.—Effect of spin rate on the deviation of the velocity vector due to a thrust misalinement angle 5 of 0.10°. 

Damped body; t~ 30 seconds (b irnout). 



Deviation of velocity vector, <r,deg 


VACUUM MOTION OF SPINNING VEHICLE 


WITH VARIABLE MASS AND INERTIA PROPERTIE: 



Kkmirk 0 Linear of thrust misalinement on the 

deviation of the velocity vector for several spin rales. 
Undamped body; / 30 seconds (burnout). 



UioeitK 10. KlTect of spin rate on the deviation of the velocity vector due to a principal-axis misalinenn 

angle e of 0.10°. Undamped body; ( 30 seconds (burnout). 





Deviation of velocity vector, cr, deg 
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ITiji'HK 


K>. -Effort of spin rati 1 on the deviation of the velocity vector due to an applied moment for a 25-iuch spher¬ 
ical rocket motor. Undamped body; / -1*0 seconds (burnout). 


CONCLUDING REMARKS 

The (‘fleets of spin stabilization on the flight- 
path angular deviation caused by initial pitching 
rate, thrust- inisalinement, and principal-axis inis- 
alinenient lmv(* heeu evaluated by an approximate, 
analytical solution which yields the flight-path 
angular deviation as a function of spin rate for a 
vehicle whoso mass and inertia properties are 
varying with time. Agreement with a numerical 
integration of tin 1 six-degree-of-freedom equations 


of motion was found to he good in all cases in 
which the inertia ratio was essentially constant. 

Spin stabilization was found to be very effective 
in averaging out the errors caused by initial pitch 
rates and thrust, misalincmcnts, but was com¬ 
pletely ineffective in controlling the errors in flight- 
path angle due to principal-axis inisalinomont. 

Langley Research ("enter, 

National Aeronautics and Space Administration, 
Langley Air Force Base, Ya., April Jl, UWi. 






APPENDIX 

EQUATIONS OF MOTION IN BODY-FIXED AXIS SYSTEM 


The symbols 
as follows: 

used in this appendix arc* defined 

FxJ'yf'z 

components of external force along 
the body AA P-, and Z-axis, 
respectively, 11) 

Ay« I r «lz 

moments of inertia about body 
A"-, Y- y and Z-axis, respectively, 
slug-ft 2 

I x z 

product of inertia, slug-ft 2 

I 

I y or J z (where I r —I x ) f slug-ft 2 


direction cosinc*s between inertial 
axes and body-fixed axes, ■! — 1, 
2, or 

U x MyM z 

external momc*nts about body A-, 
Y-, and Z-axis, rc*spc*ct ively, 
ft-lb 

m 

mass of vehicle, slugs 

]>><!,*' 

angular-velocity components ol 
vehicle about body A-, K-, and 
Z-axis, respectively, radians/see 

r„X 

longitudinal distance frojn vehicle 
center of gravity to jet exit, ft 

ir.uruu 

t ransformat ion mat ricc*s 

n } >\ir 

velocity c-omponents of vehicle 
along A-, } and Z-a\is, respec¬ 
tively. ft/SC‘(* 

V 

vehicle velocity vector, ft /sec 

XT,Z 

orthogonal body-axis system 

7,//.: 

coordinate's of center of gravity of 
vehicle along A-, ) and Z-axis, 
respectively, ft 

a lv,z 

fixed orthogonal axis system, coin¬ 
cident with body-axis system at. 
zero t ime 

7 

orientation of velocity vector in 
vortical plane, radians 

6 

pitch Euler angle, radians (except 
as noted) 

X 

orientation of velocity vector in 
horizontal plane, radians 

<T 

total orientation of velocity vector, 
radians (except as noted) 


<t> roll Euler angle, radians (except as 

noted) 

\f/ yaw Euler angle, radians (except 

as noted) 

A dot over a quantity denotes differentiation 
with respect to time. 

The rigid-body equations of motion as pro¬ 
gramed on the IBM 704 electronic computer are 
the standard six-degree-of-freedom equations of 
motion found in most aeronautical textbooks on 
dynamic, stability (for example, ref. d, eqs. (10-11)), 
with damping in pitch and yaw added to the 
second and third moment equations, respectively. 
These equations are repeated here for reference. 
The tinee force equations arc* 


m (it —rr+ w<f) - X Kx 

(A 1) 

m (r — wp + ur) — X l V 

(A2) 

m (w— ft <} + rp) — X ! 4 'z 

(A:i- 

and tin mometit equations including jet damping 
are 

/■-/)— Lx/j h (Iz —Iy) <P‘ — IxxPH — X^/a (A4 ) 

/ yH ~\ (7 v— Iz)rp-r IxAp 2 — r) 

+ l Y (i J rmri v^ Vl/r 

(A5 ) 

fyj — lx 'P 4 < / V— lx } P<l + IxxH r 

4 lyf~\ ^>\\l z 

(AOj 


The plant it ic*s used here are referred to a set of 
axes lb ed to and rolling with the body, whereas 
the am lysis in the body of this paper is referred 
to a sot of non rolling axes. 

The quantities ft, i\ w and p, <p r are found by 
integraion of these* eejuations by the* Runge- 
Kutta nethod. and are then transformed to a set 
of iner ial axes. The following assumptions are 
made* t > simplify tin* integration: 

M >ments and products of inertia vary 
in early with t ime 


IS 
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Mass varies linearly with time 

Damping coefficients (m/j x ~\-1) are constant 

The inertial quantities are found by means of 
1 1 1 e following transformation equations (ref. 4): 


The angular rates <y, r about the body axes 
are transformed to angular rates 0 , 0, 0 about the 
space axes by means of the following transforma¬ 
tions (ref. 4): 


.r 


7/ 


0 =p-ftan 0 (c/ sin 04 /* cos 0 ) 

Tf 

=r/’ii! ^ 2 ) ['Al 


CAT) 

0 — <j cos 0 —r sin 0 ^ 

0 —- —- (c/ sin 0 )-/‘ cos 0 ) 
cos 0 

- 





(Alii) 


where the transformation matrices [ 1 \\, [7\], and 
[7a] are given by 


[7i] 


cos 0 —sin 1 p 0 
sin 0 cos 0 0 

0 0 I 


As) 


r 


i 7)1 


cos 6 0 

0 1 
-sin 0 0 


l 'Al 


sin 0 

0 

COS 

0 


A9) 


0 cos 0 —sin 0 
0 sin 0 cos 0 

Multiplying the matrices gives equation (AT) as 


Alt)) 


x 

V 


/1 M I ii 1 


1 2 in 2 n 2 

U:< ":J 


U 

r 

w 


A1 1) 


where 


/. — cos 0 cos 0 
/ 2 —sin 0 cos 0 
/ :1 — — sin 0 

///,— <‘(is 0 sin 0 sin 0 —sin 0 cos 0 
/// 2 —sin 0 sin 0 sin 0 -j cos 0 cos 0 ^ (A 12 ) 
//<3 = cos sin 0 

//j = (*os 0 sin 0 cos 0-{ sin 0 sin 0 
/E» = sin 0 sin 0 cos 0—cos 0 sin 0 
{ — ms 0 cos 0 , 


Application of these transformations to the body- 
axis quantities ?/, /\ >/' and p, (/, /■ yields time 
histories of T, 7/, ” and 0 , 0. Then, from sketch 

2, 7 , X, and a can be found: 





Sketch 2. 


Thus, 


7 — tan 1 


\V'K'T 


X —tan 


a t a 11 


-1 


if 


x 
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